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Abstract 


We consider the Higgs boson decay processes and its production including all Run 
I results, through a parametrisation tailored for testing models of new physics beyond 
the Standard Model, and complementary to the one used by the LHC working groups. 
Different formalisms allow to best address different aspects of the Higgs boson physics. 
The choice of a particular parametrisation depends on a non-obvious balance of quantity 
and quality of the available experimental data, envisaged purpose for the parametrisa¬ 
tion and degree of model independence, importance of the radiative corrections, and 
scale at which new particles appear explicitly in the physical spectrum. 

The most refined constraints can only be obtained by the experimental collabora¬ 
tions at present, as information about correlation between the various uncertainties 
on the different decay modes is not completely available in the public domain. It is 
therefore important that different approaches are considered and that the most detailed 
information is made available to allow testing the different aspects of the Higgs boson 
physics and the possible hints of physics beyond the Standard Model. 



Introduction 


After Run I and with the break in data acquisition for the 2013-2014 upgrade to nominal 
energy, LHC physics experimentalists have further finalized the analysis of data that was 
gathered throughout Run I (2011-2012), which led in their preliminary steps to the announce¬ 
ment of the discovery of the 125 GeV scalar boson in July of 2012. As of April of 2015, the 
CMS and ATLAS collaboration have released most final results concerning the study at 7 
and 8 TeV of the couplings of the detected boson to the rest of the Standard Model (SM) 
particles, and are ready to start acquiring data at a higher centre-of-mass energy. 

In the Standard Model, the Higgs boson, physical particle of the Higgs doublet after the 
breaking of the electroweak symmetry, is of course coupled to the weak gauge bosons and 
fermions to which the mechanism effectively gives mass at low energies. At loop level, two 
final states experimentally stand out of the many possible Higgs decays: the one with two 
photons and the one with two gluons. Those are therefore the main interesting channels 
that are looked for in colliders, and the couplings involved can be used to constrain new 
physics models in which they are modified. Models beyond the Standard Model (BSM) 
describing the electroweak symmetry breaking often do so in a different way than the SM, 
implying a change in many of these couplings with respect to the SM predictions. Instead of 
doing computations and extracting constraints individually for each model, it is convenient 
to use model-independent parametrisations of the couplings, in order to perform the whole 
constraining procedure a few steps shorter. In (lo) , the authors advocated for the use of 
a parametrisation tailored for the study of models where loops have a sizeable impact on 
effective Higgs couplings. 


1 Higgs rescaled couplings : a loop oriented parametri¬ 
sation 

Among the particle to which the Higgs boson couples, a distinction can be made between 
particles coupling at tree level and particles couplings only through loops of other particles. 

1.1 Tree-level couplings 



g,/ R 


In the fundamental SM Lagrangian the Higgs doublet is coupled to elementary fermions 
and gauge bosons, and these couplings translate into those of the Higgs boson after symmetry 
breaking. The latter can therefore be defined straightforwardly as the coefficient in front of 
the coupling terms in the SM effective Lagrangian. 
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In a BSM model however, whether the Higgs mechanism is still relevant and additional 
physics come into play, or the whole symmetry breaking occurs in a different way, the 
Lagrangian may be different (effectively or fundamentally). Those couplings are therefore 
modified, and we can parametrise this modification through the definition of the ratio: 

9hii 
yhii 

Parameters defined in this way (in our case Kb, Kt, k^, and Ky, respectively rescaling 
the Higgs couplings to down-type quarks, up-type quarks, charged leptons, and weak vector 
bosons) obviously have a value of 1 in the SM, and may deviate from this value when the 
model is different from the SM. 


1.2 Loop-induced couplings 



In order to use the experimental constraints from those channels to test a model, it is also 
convenient to use simple parameters allowing comparison between experimental results and 
model predictions. However, as stated previously, the Higgs couplings to pairs of photons 
or of gluons do not appear explicitly in the fundamental Lagrangian but only appear in an 
effective description when considering the effects of loops. Therefore, the parameters are 
trickier to define in a unique way. One of the standard approaches is to do the same as 
in the case of direct couplings and rescale the effective couplings of the model to the same 
effective couplings in the SM. 

However, this approach does not distinguish the variation of this effective coupling due 
to the change of the direct couplings from the one due to additional loops coming from 
the new physics. Another approach, introduced in is interesting in the case of models 
where additional particles appear. In this case, one can take into account the loop-induced 
character of these couplings and introduce the change in the effective coupling at the level 
of the loop calculation. More precisely, in the SM one can compute the partial width of the 


Aw{tw) + 3 At{Tt) + ... 

+ ... 


Higgs in those channels through; 
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where A^/ and At represent the analytical amplitude of the W boson and top quark 
loops, Ct and Cf include QCD corrections, and the dots represent negligible lighter fermions 
contributions. Dehning r = the relevant expressions are: 


M^) 

= ^{T+{T-l)f{T)) , 

Aw{r) 

= -;^( 2 r 2 + 3 r + 3 ( 2 r-l)/(r)) , 


( arcsin^^r t < 1 

and /(r) 

I i [log ^ > 1 


Additional loop diagrams coming from new physics contribute to these amplitude, and 
we therefore dehne the parameters and Kgg as: 
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which gives for the partial widths in a new physics model: 

Gpa^m\ 

"" ^ 128v^7r3 
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16v^7r3 

The hrst thing one can notice here is that the contributions to the variation of the width 
from the modihcation of the tree-level coupling and from the existence of new loops are here 
decorrelated, since the hrst are taken into account thanks to the Ky and nt parameters. On 
the other side, the normalisation to the top amplitude is arbitrary, but stays general and 
might gives interesting values for models where the new physics include top partners. 

As a side comment, we can note that in the case where the custodial character of the 
electroweak symmetry breaking is questioned (i.e. when couplings to the W and Z bosons 
are no longer rescaled by a common parameter), there is an advantage in considering the 
separate amplitude calculation rather than a simple width rescaling, as the parameter rescal¬ 
ing the amplitude of the W loop can then be described as kw In this prospect, additional 
information can be obtained in this channel from the use of this parametrisation, as kw 
appears explicitly in the expression of T^.^, as opposed to kz- However, as we will see in 
the second part, most LHC published results use a common rescaling factor for W and Z 
bosons couplings, making this information impossible to recast here, as it is only accessible 
to experimentalists. 
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On another hand, the parameter describing the loops coupling the Higgs boson to a pair 
of gluons can also be used to describe the modihcation of the cross section of production of 
Higgs through the fusion of gluons from the colliding protons, which is the main production 
mode at hadron colliders. We consider here that the effective coupling is the same for 
production cross section and decay width, although this might not be the case due to PDF 
effects. 

1.3 Interpretation of BSM models 

In particle physics experiments, measurements rely on the counting of events, therefore 
depending on differential cross sections, which can often be calculated in given models. As 
a consequence, the couplings involved in the calculation of cross sections have an influence 
on the expected result. 

More specihcally, one can dehne the signal strength /i as the number of events expected 
in a model for a given signal normalized to the number of events expected in the SM; 

= jj, . They are therefore functions that depend on as many parameters as the 
cross sections. However, if the additional physics is expected not to change the shape of 
the cross-section distributions, it is reasonable to consider only one constant signal strength 
per measured physical process. Same as the cross sections, the signal strengths depend on 
the values of the couplings, and we can therefore write them as functions of our rescaling 
parameters: 


NP I \ NP 

, ^ fl yKv J Kt, Kh, Kq-J K-fy J Kgg) O' 

/r [Kv , Kt, Kb, Kr, Kyy, Kgg) = “ ^SM 

2 Recasting Higgs constraints 

2.1 LHC data 

In the case of the Higgs boson analyses, the previous dehnitions would provide us with one 
signal strength per combination of production mode and measured decay channel. 

However, several assumptions make this number decrease. First, under the narrow width 
approximation (which experimentally seems valid), the cross sections for a given process can 
be factorized as the Higgs production cross section times the branching ratio of the hnal 
state. We can therefore use only one rescaling factor per production mode and one per 
decay channel. Furthermore, as VBF and VH production modes rely at tree level on the 
Higgs coupling to weak gauge bosons, we can expect them to be rescaled in the same way. 
On the same level of idea, in the SM the ggh production mode relies mostly on a top quark 
loop, and therefore scales in the same way as the tih production mode. Although this is 
not exactly true (or we wouldn’t have included a Kgg parameter), there is no analysis at the 
moment that include both the ggh and tih production modes (the tih production mode is 
only probed through Higgs decays to hh). 
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One of the forms under which the experimental collaborations provide the results of 
their analyses is as exclusion contours for given conhdence levels, in the plane mapped by 
the production modes signal strengths Hvbf/vh and Hggh/tth- Each studied hnal state gives 
an exclusion contour, corresponding to the constraints on the signal strength of the whole 
process, from production to decay. 

Fig. 5 in 0 summarises the hnal results of the different analyses for 7 and 8 TeV of 
the CMS collaboration, gathering the various exclusion contours at 68 % C.L. for the various 
decay channels. 

In the case of the ATLAS collaboration, a single hgure gathering the contours obtained 
from each analysis after the legacy analysis of Run I was not published (as of writing of this 
paper). Therefore, the individual relevant hgures were gathered from the analysis papers 
( 77 , Fig. 20 in j^, Fig. 12 in [^, WW* Fig. 40 in |5) and ZZ* —)■ 4£ Fig. 20 in [^) 

As an addition, ATLAS results from the bb Higgs decay searches were included as an 
exclusion band around a central value (Fig. 20 in Q). This comes from the fact that this 
analysis relies on a selection tailored for events where the Higgs boson is produced through 
the VH production mode, using the presence of the additional gauge boson to reduce QCD 
noise. As a consequence, this analysis does not constrain HggH/ttH- 

2.2 Fitting contours as ellipses 

As these exclusions take into account possible statistical fluctuations, experimental and the¬ 
oretical systematic uncertainties on the measurements, one can assume that the likelihoods 
of the signal strengths, for large numbers of measured events, are distributed as Gaussians. 
This means that the log-likelihood functions are paraboloids. In terms of constraints, this 
can be translated into the fact that the exclusion contours at a given conhdence level are 
ellipses. 

We therefore htted the parameters of each log-likelihood paraboloid, so that the positions 
of the points of the contours taken from the plots matched those of the ellipsis of points with 
hxed log-likelihood equal to 2.3 , i.e. the value at which lays the 68 % C.L. of a 2 degrees of 
freedom (d.o.f.) distribution. The comparison of the LHC results and their htted ellipes 
is presented on hgure In order to further test the validity of the hypothesis of gaussianity 
of the log-likelihood functions, we compared the 95% C.L. exclusion contours extrapolated 
from the htted distributions to the experimental contours from the ATLAS experiment. 
The comparison plot is displayed in hgure and show reasonable agreement. 

2.3 Constraining models 

After obtaining the shape of the likelihood for the signal-strengths, they can be translated 
into a single likelihood for the parameters dehned previously. However, the proper combi¬ 
nation of the information from the diherent channels can only be done by experimentalists, 
as systematic uncertainties might be correlated, through theoretical uncertainties in the de¬ 
scription of the same decay channel or production mode, even between diherent experiments, 
and through experimental systematics, within a same experiment. Being aware of this, we 
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Figure 1: Relevance of the hypothesis of a Gaussian likelihood: comparison of experimental 
68%C.L. exclusion contours with the fitted contours (right; CMS data, left; ATLAS data). 
Colour code per final state: purple: bb, yellow; 77 , green: rr, red: WW, blue: ZZ. Gray 
rhombus; SM 


nevertheless combine the different channels by multiplying the various likelihood functions. 
Moreover, as the experimental contours can only be interpreted as likelihood-ratios, the most 
likely point being the null hypothesis, we combine our different functions as such, giving us 
the likelihood ratio function C, or alternatively the log-likelihood ratio function Ay^: 


C {Kj) 


or (nj) 


rii i^j) 

rii (%) 


where Kj represents the set of all k’s, and kj their values maximising the total likelihood 
(minimising the total log-likelihood). 

According to Wilks’s theorem, the log-likelihood function follows a y^ distribution with 
a number of d.o.f. equal to the number of parameters on which the function depends. In our 
case, with 6 parameters, the Ay^ function should therefore follow a 6 d.o.f y^ distribution. 

In order to see how this applies to specific models, we provide a few examples of models 
that can be described by our parametrisation, and for which values of the parameters have 
been described in [l 0 | : 

■ Colour octet model (^: Model where the scalar SM sector is extended with a colour 
octet. As this octet can be decomposed into SU(2) representations, additional neutral 
and charged scalars appear. The model point represented in Fig. corresponds to 
ms = 750, Ai =4, A 2 = 1, and A 3 = 2 A 2 , while the line corresponds to a varying 
value of ms- 
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Figure 2: Relevance of the hypothesis of a Gaussian likelihood: comparison of experimental 
95%C.L. exclusion contours with contours extrapolated from the 68%C.L. ht (ATLAS only). 
Colour code per hnal state: yellow: 77 , green: rr, red: WW, blue: ZZ. Gray rhombus: SM 



_ 3 Xiv^ ^(8) (2Ai + 

~ 2 ^^ '""" “ ^ R4 

(g) 5D UED l^: A Universal Extra Dimension model with one extra dimension, where 
only the top and W Kaluza-Klein resonances contribute and the result scales with the 
size of the extra dimension. The model point here corresponds to rriKK = 500 GeV, 
and the line corresponds to a varying ttikk- 


63n^ / mw \ ^ / rut 

16 X 6 \mKK) 6 \mKK 


^aa 


^ / mt y 
6 Xtyikk) 


A Simplest Little Higgs |^: The electroweak SU{2) gauge group is in this model em¬ 
bedded in a larger SU{3), spontaneously broken by triplets through a Higgs-like mech¬ 
anism. The corrections scale with the W mass, and the point corresponds to the value 
mw' = 500 GeV within electroweak precision constraints for a model with T-parity. 
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* Littlest Higgs [^: The Higgs in this model is one of a set of pseudo-Goldstone bosons in 
an SU (5)/S'0(5) nonlinear sigma model. The result scales with the symmetry breaking 
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scale /, which can also be set to low values for a model with T-parity (there is 
also a mild dependence on the mass of the extra gauge boson contributing to 
arbitrarilly hxed to be equal to / and on the triplet VEV x, that we set to x = 0); 
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Figure 3: Constraints on various models, 
All other parameters are taken equal to 1. 
68%, 95% and 99.7% C.L. values for a 6 d.o.f. 
suits from former LHC analyses, and the right o 



Kyy 

in the plane spanned by n^-y and Kgg. 
The distribution is compared to 

distribution. The left hgure uses re- 
ne uses results from the legacy analyses. 


Figure represents the contours of the regions excluded by two sets of LHC results, at 
68%, 95% and 99.7% C.L., in the plane spanned by Kgg and The SM is represented by 
a grey rhombus, at the centre of the hgure. The right plot was produced using legacy data 
described before, whereas the left one was produced using the same LHC results as in and 
is displayed here to see the evolution of the constraints. We can see that the latest bounds 
are more stringent but also shifted so that the centre is closer to the SM. Therefore, some 
values of the parameters that were excluded at a given conhdence level no longer are, as we 
will see a bit further. 

More precisely, from these constraints, we can compute values of the parameters for which 
the models are excluded, for instance along the lines drawn on the hgure. In table four 
values are provided, corresponding to two diherent approaches of the way the model should 
match the observations, each one both for former constraints and the recent legacy ones. 

The hrst and third lines of the table show the values of parameters excluded when com¬ 
paring the model’s likelihood Ay^ to a 6-dof y^ distribution. This means that we consider 
that the model in itself does not describe all new physics, and other modihcations might 





















have an influence on Hb, •••, but we slice the parameter space to observe the plane where 
these parameters are 1. We then get a parameter value under which the model is excluded 
at 95% C.L.. 

On the second and fourth lines, we instead consider that the only modification to the 
Higgs couplings come from the model we are testing. The other k’s are then no longer free 
parameters, and Wilks’s theorem asserts that the log-likelihood ratio function now follows a 

distribution with n d.o.f., n being the number of parameters of the model (and the null 
hypothesis being of course the total minimum, amongst points in k space that can be 
described by model). The bounds are therefore tighter, and the values excluded higher. 




Colour Octet 

5D UED 

Simplest Little Higgs 

Littlest Higgs 

Previous 

analyses 

6 d.o.f 

ms = 1040 GeV 

fTT-KK = 480 GeV 

mw' = 200 GeV 

mw' = 475 GeV 

n d.o.f 

ms = 1220 GeV 

'tokk = 660 GeV 

mw' = 280 GeV 

mw' = 540 GeV 

Legacy 

analyses 

6 d.o.f 

ms = 835 GeV 

m-KK = 430 GeV 

mw' = 250 GeV 

mw' = 580 GeV 

n d.o.f 

ms = 925 GeV 

'n^KK = 530 GeV 

mw = 380 GeV 

my[/i = 690 GeV 


Table 1: Values of the parameters excluded for different models under specific assumptions. 
The first two lines correspond to old constraints and the two last lines to legacy constraints. 
In each case, in the first lines model predictions are interpreted as following a 6-d.o.f 
distribution and in the second line as a y^ distribution with number of parameters equal to 
the number of parameters of the model 


From the values in this table, we can note that the most stringent bounds for given 
models are not necessarily the most recent ones. For instance, the colour octet and the 
UED models are less constrained by legacy results than by the former constraints. This can 
be easily understood by comparing the two diagrams in figure The constraints on Little 
Higgs models, on the other hand, follow the expected variation, and are more constrained 
by more recent results. 


Conclusion 

The parametrisation of the Higgs boson couplings described here, laying between the simple 
parametrisations with few parameters and the more complete but complex ones, is com¬ 
plementary to other similar parametrisations, and could be used in its full potential if 
constrained directly by experiments. In the meantime, the procedure presented here is a 
straightforward way for theoreticians to get constraints on models including new loops in 
Higgs physics. As an example, we computed bounds on a set of models with constraints 
extracted from CMS and ATLAS Run 1 Legacy analyses, showing this parametrisation is a 
simple and effective way to constrain BSM models. 
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Note: The ATLAS collaboration released during the writing of this paper a conference 
note containing a figure summarising all Higgs signal strengths constraints. It can be found 
in Fig. 3 in Q, however, as this figure is marked as a preliminary result, and we did not 
expect the results to be different from the final analysis papers, we did note include it in our 
study. 
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